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Abstract. In this paper, the nonlocal mixing of coherent current states in d-wave 
superconducting banks is investigated. The superconducting banks are connected via 
a ballistic point contact. The banks have a mis-orientation and a phase difference. 
Furthermore, they are subjected to a tangential transport current along the a6-plane 
of d-wave crystals and parallel with interface between superconductors. The effects 
of mis-orientation and external transport current on the current-phase relations and 
current distributions are subjects of this paper. It observed that, at values of phase 
difference close to 0, tt and 2tt the current distribution may have a vortex-like form in 
vicinity of the point contact. The current distribution of this junction between d-wave 
superconductors is totally different from the junction between s-wave superconductors. 
As the interesting results, spontaneous and Josephson currents are observed for the 
case of 4> = 0. 
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1. Introduction 



The weak link between <i-wave superconductors is a long studied problem theoretically 
[HI21E1H1IS1IE10IE1I3]- A theoretical investigation of the total transparent Josephson 
junction between two (i-wave superconductors has been done in Ref. p. Using the 
quasiclassical approach, a (i-wave Josephson Junction with low-transparent interface 
has been studied in Ref . p] . Anisotropic and unconventional pairing symmetry has been 
considered for d — / — d systems and zero energy states (ZES) as the result of sign 
change of the order parameter have been observed in paper [3]. A spontaneous current 
parallel to the interface between <i-wave superconductors has been found in paper pQ. 
The junction between current-carrying states of d-wave superconductors, has been 
investigated in Ref. pi]. Authors of paper pE] by numerical self-consistent calculations 
shows that, the supercurrent parallel to the junction may flow in the direction opposite 
to the current direction at the superconducting banks. Effect of transparency of interface 
of (i-wave Josephson junction was studied in Ref. [5]. In Ref. [6], effects of transparency 
and mis-orientation of two d— wave crystals have been investigated analytically. Using 
the Bogoliobov-De Gennes equations ZES as the origin of zero bias conductance 
peak (ZBCP) were studied in Refs.0 Ej. ZES are introduced as the fingerprint of 
unconventional pairing symmetry in[71 [8]. In Ref. [9], a special geometry of d— wave 
superconducting layer as a weak link has been investigated and the n Josephson junction 
has been observed. Also because of high critical temperature of (i-wave superconductors 
(and cheap production technology), many experimental works about the <i-wave weak 
link have been done in the las two decades pm [HI [121 EH O 13 [T7]. A complete 
review of these experiments has been presented in a review paper [10J. A phase 
sensitive experiment (phase of the superconducting order parameter) has been presented 
by authors of [11] for determination of the symmetry of order parameter in high T c 
cuperate superconductors. Using the phase interference experiments in the Josephson 
junctions, d— wave pairing symmetry in the cuperate superconductors has been observed 
in Refs. [TT1 [T2| [T3] . A nonsinusoidal form of current-phase diagram has been observed in 
[H] , experimentally. Authors of Ref. [15] have measured the current-phase relationship 
of symmetric grain boundary weak link and observed that, temperature controlled sign 
change of the first harmonic of the Josephson current (I(<p) = I\ sin <fi + I2 sin 2<p + ■ ■ ■). 
Because of competition between the first and second harmonics of Josephson current, a 
nonmonotonic temperature dependence of the critical current has been reported in [15J. 
An experimental investigation of a Josephson junction between d— wave superconductors 
has been done and effect of insulator between them has been studied in [17J. They 
observed — 7T transitions by reducing the width of insulator in d — I — d Josephson 
junction. From the other hand, it is well-known that, non-locality and Josephson effect 
are coexisting. Charged particles orbiting around a magnetic flux are influenced by 
a magnetic flux as a phase difference although, the region including flux is forbidden 
for charged particles. This phase as the Aharonov-Bohm phase is a demonstration of 
non-locality of quantum mechanics. While the supercurrent in a superconducting bulk 



Effect of Parallel Transport Currents on the D-wave Josephson junction 



3 



depends on the phase gradient locallyJTS], j(r) oc Vy(r), the Josephson supercurrent 
depends on phase difference non- locally [19], j(<^ 2 — y?i) oc sin ((p 2 — (fx) oc ((p 2 — <fi). 
The interplay between local supercurrent states and non-local phase difference between 
the superconducting bulks is an interesting problem. Local supercurrent states are 
introduced by superfiuid velocity of Cooper pairs. The, non-locality of Josephson 
current in the point contacts and the effect of super-fluid velocity on the current 
states in narrow films and wires have been studied in Ref. [20j. An anomalous periodic 
behavior in terms of magnetic flux has been observed in Ref.|20j. This anomalous 
property is demonstrated as a result of a non-locality of supercurrent in the Josephson 
junction[20j. In addition, experimental results of Ref. [20] have been confirmed in 
analytical calculations of Ref. [21]. The dynamical Josephson junction between s- 
wave superconductors has been investigated in Ref. [22J. They studied the quantum 
interference between right and left s-wave superconductors, in which parallel transport 
currents are flowed. The existence of two anti-symmetric vortex-like currents near the 
contact and at <fi ~ ir as a new phenomenon was reported in Ref. [22]. The authors of 
Ref. [22J found that the total current is not the vector sum of Josephson and the transport 
currents because of a new term in the current. This term is called the "interference" 
current and can be a "parallel Josephson current". In the Ref. [22], effect of reflection at 
the interface between s— wave superconductors has been investigated analytically and 
numerically. 

In this paper, a planar weak link between two d-wave superconductors with a 
phase difference between order parameters is investigated. The afr-planes of two 
superconductors have a mis-orientation and c axis of two crystals are parallel to the 
interface between d-wave superconductors. In the center of the interface we create 
an ideal transparent thin slit with length L and width a. Interference between wave 
functions of the left and right superconductors occurs through the slit. The remainder 
part of interface is an ideal insulator and is impenetrable for Cooper pairs. Also, two 
transport supercurrents at the afr-planes are flowing parallel to the insulator and the 
contact plane (see Fig{T]). The Josephson current from one of the bulks to another, is 
a result of the interference between states with phase difference as is predicted in 
Ref.[19j. The contact scales (thin slit), length L and width a, are much larger than the 
Fermi wavelength and smaller than coherence length of superconductivity. Furthermore, 
these scales are small as compared with the mean free path of quasi-particles. Therefore 
the quasi-classical approximation for the ballistic point contact can be used. The 
Eilenberger equations for this system are solved and Green functions are obtained. The 
effects of mis-orientation and phase difference between order parameters and super-fluid 
velocity on the current distributions and current-phase graphs are investigated in this 
paper. 

The organization of the rest of this paper is as follows. In Secj2] the quasiclassical 
equations for Green functions are presented. The obtained formulas for the Green 
functions are used to analyze a current state in the ballistic point contact. Also the 
effects of transport current and mis-orientation on the current distribution at the contact 
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Figure 1. Model of the contact in the insulating partition, along the y, between 
two mis-oriented d-wave superconducting bulks with transport supercurrent on 
the banks. The plane of paper is afe-plane of d-wave superconductors. In the 
d-wave superconductors like YBaCuO the a6-plane is the plane of CuO. 

plane are investigated. In Secj3] the results of simulation for current distribution in 
vicinity of the contact will be investigated. An analytical investigation of system near 
the critical temperature will be done in SecJU The paper will be finished with some 
conclusions in Secj5j 



2. Formalism and Basic Equations 

The Eilenberger equations for the ^-integrated Green's functions are used to describe 
the coherent current states in a superconducting ballistic micro-structures [23J: 

d ~ ~ 
vf • ttG,(v f , r) + \uj% + A(v F , r), G u {v F , r)] = (1) 
or 



where 
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Pf V T c = 0.0 
Pf V T c = 0.22 
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Figure 2. Tangential current j y versus <j) for T/T c = 0.1, a\ = a r = 

in the units of jo = 4:ireN(0)vFT. This is like the case of junction between 
conventional superconductors [22] 



A is the superconducting order parameter, f 3 is the Pauli matrix, and (^(vp,! - ) is the 
matrix Green function which depends on the electron velocity on the Fermi surface vp, 
the coordinate r and Matsubara frequency uj = (2n + l)vrT, with n and T being an 
integer number and temperature respectively. Also the normalization condition 

g u = yfi - Ut (3) 

with fl being time-reversal counterpart of f w should be satisfied by solutions of the 
Eilenberger equations. In general, A depends on the direction of vp and r and it can 
be determined by the self-consistent equation 

A(v>, r) = 2vriV(0)T ]T (V(v F , V F )/ u (v' Fl r)) / (4) 

LU>0 

and the current density by 

j(r) = -47rieN(Q)Tj2(vF9Uv F ,v)) VF (5) 

respectively, where V(vp,v f F ) is the interaction potential, iV(0) is the 2D density of 
states at the Fermi surface for each spin projection and (...) is the averaging over 
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Figure 3. Tangential current j y versus phase 4> for T/T c = 0.1, ai = and 

a r = -in- 
directions of vp. Solution of the matrix equation together with self-consistency 
gap equation (j4]), Maxwell equation for superfluid velocity and normalization condition 
determines the current j(r) in the system. The thickness of c?-wave superconductors are 
assumed to be smaller than coherence length £ = • Thus the spatial distributions of 
A(r) and j(r) depend only on the coordinates in the plane of the film and Eilenberger 
equations (pQ) reduce to 2D equations. 

Also equations ([I]) for Green functions G w (r,vp) have to be supplemented by the 
condition of specular reflection at the region (x = 0, \y\ > a) and continuity of solutions 
across the point contact (x = 0, \y\ < a). Far from the contact the Green functions 
should be coincident with the bulk solutions and the current should be homogeneous 
transport current along the y-axis. 

In this formalism, the current has to be determined by a self-consistent gap equation 
fll]) together with the Maxwell equation for superfluid velocity (Ampere law d ^f 3 ^ = 
—fioJy 0t (x) with jiQ for free space and v s = —^A y (x) for superfluid velocity as in paper 
[5J [21] ). While for simplicity, the self-consistency of order parameter is ignored and a 
step function is considered for spatial dependence and we do not consider effect of current 
distribution on the superfluid velocity. We believe that as in the papers [251 126], the 
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self-consistent investigation of d-wave Josephson junction does not show a qualitative 
different from the non-self-consistent results. 

For A and v s being constants at each half plane an analytical solution for Eilenberger 
equations can be found by the method of integration along the quasi-classical trajectories 
of quasi-particles. In any point, r = (x,y), all ballistic trajectories can be categorized 
as transit and non-transit trajectories. 

For the transit trajectories the Green functions satisfy continuity at the contact and 
the non-transit trajectories satisfy the specular reflection condition at the partition, 
(x = 0,\y\ > a). Also, all transit and non-transit trajectories should satisfy the 
boundary conditions in the left and right bulks. Making use of solutions of the 
Eilenberger equations, the following expression is obtained for current at the slit: 

■i ni 1/ a \ a Mfn\ ^V" /- tijtti + Qr) -ir/A / A 7 .sin0 \ 

j(x = 0, \y\ < a, 0, v s , a h a r ) = 4ireN(0)v F T ^ I v — — ) (6) 

~\ iliil r + a? + A; A r cos <fi 

where, A/ /r = Ao cos(2(0 — cx^ r )) for d x 2_ y 2 symmetry Vt^ r = ^Ju 2 + Af >r , uj = oj + ippV s 
with uj being Matsubara frequency and v s is super-fluid velocity, v = vf/vf is the unit 
vector, r] = sign(v x ). In this non-stationary Josephson junction, v s ^ 0, the current 
has both } x and } y components. We define the Josephson current, external transport 
current, spontaneous current and interference current as: 

iJosephson = j(0, V a , a r )x (7) 
'^Transport = j (0 = 0, V s ^ 0, a t = CX r = 0) y = '^Bulk, (8) 
^Spontaneous = j (4> °, V s = 0, «/ ^ 0, a r ^ 0) y (9) 
j JT — j(0> v si tt/, a r )^ — j Spontaneous — j Transport (10) 

respectively. The Josephson current, jj = j x , is normal to the interface between 
superconductors as was considered by B. D. Josephson in [19] and parallel component of 
current j y consists of three terms of current, an external transport current, spontaneous 
current and "interference" current. The new current term, "interference" current, 
depends on the super-fluid velocity, orientations with respect to the interface and phase 
difference between order parameters. This term of current, jjr, is completely different 
from the transport current on the banks jr. 

Thus in addition to the spontaneous current for the stationary <i-wave Josephson 
junction that is investigated in Ref. [6] and transport current, we observe another current 
parallel to the interface (jjr)- In particular, at (f> ~ n it may go into the opposite 
direction to the external transport current on the banks (depending on the orientations). 
This sign reversal of tangential supercurrent which is origin of vortex-like currents near 
the orifice has been observed already in Ref.[l] for the case of d— wave junction. At 
two sides of vortex-like currents, currents are flowing parallel and antiparallel to the 
external supercurrent at the superconducting banks. In the Ref.jl] a superconductor- 
normal-superconductor trajectory for particles has been considered while there is one 
superconductor coupled the normal metal. So, because of one superconductor in the 
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structure of Ref. [4J it was impossible to consider the phase difference. In the case of 
junction between s— wave superconductors in Ref.[22j, it was observed that sign reversal 
can be seen only for ~ tt. So, for s— wave counterpart of setup of paper [I] it is 
impossible to see the sign reversal because phase difference does not meaning for 
system including one superconductor. While in our calculations in this paper we have 
observed that, this sign reversal and vortex appearance can be seen at (f> = 0, depending 
on mis-orientation. For the case of d— wave in paper of [I] in the absence of phase 
difference, a sign reversal of current has been observed as well as our results for suitable 
mis-orientation of ab— planes(a^ = 0, a r = | or planes (1,0,0), (1,1,0)). The planes ai = 
and a r — | are corresponding to planes (100) and (111) respectively. 



Pf V T c=o o 

Pf V T c =0.22 




Figure 4. Tangential current j y versus phase 4> for T/T c = 0.1, a\ = and 
a r = \. 

3. discussion 

The effects of mis-orientation and super-fluid velocity on the current-phase relation and 
current distribution are investigated numerically and the results are following. 
1) The vortex-like currents appear at values of <fi when the parallel current is maximum 
and positive while the external transport current is negative (Figs J2ll3lHl and 151) . 
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PfV t =0.0 
PfV T =0.22 




Figure 5. Tangential current j y versus phase for T/T c = 0.1, ai = and 

Oir = "f . 

2) At «| = a r = and arbitrary 0, the current distributions and current-phase 
graphs are identical with the s-wave which is investigated in Rcf. [22\. For example 
at ct\ = a r = and = 7r two anti-symmetric vortex-like currents are observed and 
their common axis is normal to the interface (Figj6j). 

3) At on = 0, a r = ±| and = n, two vortex-like currents are observed. Their common 
axis is rotated as much as =F~ (FiglHland Figj9]). This case occurs for a\ = 0, <p — 71 an d 
arbitrary a r with the rotated axis rotating as much as (— a r ). 

4) The appearance of vortex-like currents can be controlled by mis-orientation. For 
example, at = and a\ = a r = 0, we cannot observe the vortex-like currents as we 
cannot observe in s-wave junction (Fig|7]). However, for = 0, a^ r = 0, a r j = | the 
vortex-like currents appear since from ([6]) and for d x 2_ y 2 symmetry we have : 

7T 

j(r, = 7r, v s , a h a r ) = j(r, = 0, v s , a h a r + -). (11) 

5) In Figs l3f4f5l it is observed that for = and consequently without any external mag- 
netic flux, the interference between coherent current-states can occur and the vortex-like 
currents can be observed because mis-orientation plays the role of the phase difference 
in Eq.tfTTi). 
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Figure 6. Vector plot of the current for T/T c = 0.1, Ppv s / Aq(0) = 0.5, = 7r, 
ai = a r = 0. Axes are marked in the units of a. Because of our non-self- 
consistent formalism £o = 5a. It is similar to the s-wave junction in [22J. 

6) The parallel current, j y , is plotted in terms of the phase difference for different super- 
fluid velocities and at ct\ = a r = 0, the maximum value of current and appearance of 
vortex-like currents occur at = tt. In this case, far from = tt, we observe a constant 
current that is the external transport current on the banks (Figj2]). 

7) For a^ r = 0, a r j — ±| the maximum values of the parallel current, j y , and conse- 
quently the vortices appear at = 0, = tt, = 2tt (Figs j4ll5|) . 

8) For ai tr = 0, a r> i = | the current-phase graphs are similar to the a>i = a r = but 
a displacement as much as tt occurs (Figs J2ll3l) . Thus the vortices can be observed at 
= and = 2tt. But at = tt, we do not observe the vortex-like currents. This can 
be another difference between conventional and unconventional Josephson junctions. 

9) The superposition of dash lines in Figs l4ll5l for zero super-fluid velocity, and all the 
lines of Figj2j for zero mis-orientations gives us the three lines of FigsHj5] apparently. 
This means that in this case the tangential current is algebraic sum of transport current, 
"interference" current [22] and spontaneous current [Sj. 

10) The tangential current for a; = a r = and a^ r = 0, a r j = f is an even function 
of Josephson phase but for a^ r = 0, a r j = ±| it is neither even nor odd function of 
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-10 1 

Figure 7. Vector plot of the current for <j) = 0, a.\ = 0, a r = and T/T c = 0.1, 
Pfv s / Aq(0) = 0.5. Vortices disappeared but transport supercurrent flows. 

phase difference and the symmetry will be broken. 

The superfluid of pairs creates the transport current but the spontaneous current is 
produced by mis-orientation plus to the phase difference. However the "interference" 
current depends on all of the parameters (phase difference, mis-orientation and super- 
fluid velocity) and it can be the result of the non- locality of supercurrent. 
The spatial distributions of the order parameter and the current near and precisely 
at the contact are calculated using the Green function along the transit and non- 
transit trajectories numerically. Transit trajectories for each point are coming from 
the orifice (transparent part of interface x — 0, \y\ < a) while non-transit trajectories 
form the remainder part of interface which is impenetrable (reflective part of interface 
x = 0, \y\ > a). The current distributions are calculated and simulated numerically 
for T/T c = 0.1, pfv s /Ao(0) = 0.5 and different choices of mis-orientation for <fi = n 
(Figs IB1l6ll9l and = (FigJTJ). Following the Ref. [22], the "interference" current that 
plays a central role for production of vortex-like current is found. At the phase values 
O<0<|, ^tt < (p < 2n and ai = a r = 0, the interference current is very small thus 
the total current is close to the vector sum of transport, spontaneous and Josephson 
currents. But, for | < < |7r the "interference" current appears and the total current 
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Figure 8. Vector plot of the current for <j) = tt, ai = 0, a r = j and T/T c = 0.1, 
Pfu s /Ao(0) = 0.5. Axes are marked in the units of a. 

deviates from the vector sum of the three old currents [22\. Also, for zero mis-orientation 
the spontaneous current is zero. 

The "interference" current is always anti-parallel to the transport current. But the 
spontaneous current may be parallel or anti-parallel to the transport current. If alge- 
braic sum of transport, "interference" and spontaneous current is anti-parallel to the 
transport current on the banks, we can find the vortex- like currents (Figj6j). 
Thus the appearance of the vortex-like currents can be controlled by mis-orientation 
and phase difference. It is remarkable that, far from the contact |r| ~ £o > Q for all 0s, 
mis-orientations, temperatures T and superfluid velocities, the distributions of currents 
tend to the tangential transport currents on the banks. 

4. near the critical temperature 

For temperatures close to the critical temperature, T C — T^T C problem is solvable pure 
analytically [22]. At the contact we have 

j = jj + ispont + k + jjT (12) 
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Figure 9. Vector plot of the current for eft = ir, o.\ = 0, a r = — ? and T/T c = 0.1, 
Ppu s /Ao(0) = 0.5. The axis of vortices is rotated. 



jj = 2j c sm<f)(v x sign(v x )(^^-)j (13) 
jspont = 2j c sin0( v y sz#n(^)(^^)\ (14) 



-A/A 



45 , v 



jjT = j c A;(l - cos0) (v{;j / (^^)\ (16) 



A^ 



V 



where, as in Ref. [22] 

2 



jt(riiij) ^VM (17) 



is a critical current of the contact at (T c — T) T c , is a standard notation 

fc = (14^(3)/7r 3 )( W /T c ). (18) 
For the high value of temperatures near the T c , the critical values of currents have a 



linear dependence on the A , which can be replaced by A = \j ( 21^3) j T C (T C — T). 
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On the other hand the spontaneous and Josephson currents are sinusoidal function 
of phase difference as is expected for the currents near the T c , in spite of these two 
terms of current the "interference" current is an even function of the phase difference. 
The current is divided into the four parts, Josephson current jj, spontaneous current, 
transport current in the banks jy and the "interference" current jjr- It is observed that, 
the currents generally and near the T c obviously, depend not only on the mis-orientation 
| a; — a r \ but also depend on the orientations with respect to the interface. Because in 
the expressions (vA;A r ...), the result of angular integrations may include both \ai — a r \ 
and | a; + a r \ terms. For example, for Josephson and spontaneous currents in (fT3|) and 
(|T4"1) by angular integration on the Fermi surface we have 



respectively. Also it is obtained that for = tc and exactly at the contact, the "inter- 
ference" current jjr is anti-parallel to the jr- For = tc the "interference", Josephson 
and spontaneous currents are jj^ = — 2jr, jj = and j spon t = respectively, conse- 
quently it is obtained that j y = jr +jjT+h P ont = —Jt- In this case while the Josephson 
current is zero the terms j y and ]t that are directed opposite to the each other, control 
the appearance of vortex-like currents in vicinity of the point contact. In addition, for 
= 0, ai >r = and a r> i = |, we can observe the vortex-like currents. This property can 
be a difference between <i-wave and s-wave Josephson junctions. Because in the s-wave 
Josephson junction the vortex-like currents are observed only for = tc (Ref.[22j), while 
in the present calculations for d-wave Josephson junction, the vortex-like currents may 
appear even for = 0. Thus, for fixed values of the temperature and superfluid veloc- 
ity, the presence of vortex-like currents can be controlled by mis-orientation and phase 
difference. 



5. Conclusions 

The vortex-like currents in vicinity of the point contact are observed for c?-wave 
Josephson junction as well as for s-wave junction in Ref.[22]. The interference current as 
a result of non-local supercurrent states, appears. It may flow opposite to the external 
transport current. For = tc, ai = and a r = ±| the vortex-like currents with 
the rotated axis are observed. But as is obtained in Ref. [22], the axis of vortex-like 
currents in the s-wave Josephson junction is normal to the interface. Thus, this rotated 
axis can be used to distinguish between s-wave and d-wave junction. Also it can be 
exerted to distinguish between the junction between two pure s-wave and mixing of 
conventional and unconventional order parameters (eg.d + is). In addition to the ZES, 
this rotated axis can be another "fingerprint" of <i-wave pairing symmetry. Another 




(19) 



and for spontaneous current 

, —Sj c sin 




(20) 
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interesting result is, the behavior of system in the absence of the phase difference. For 
s-wave system only in the presence of the phase difference, <fi ~ tt ^ 0, the vortex-like 
currents appear [22], while for <i-wave Josephson junction at zero phase and zero external 
magnetic flux, it is possible to observe the vortex-like currents for some mis-orientations. 
This can be a theoretical reason that, the mis-orientation plays role instead of the phase 
difference (Josephson phase). In the stationary Josephson junction v s = the tangential 
interference current (spontaneous current) will be observed only for <fi ^ but in the 
opposite case v s ^ the tangential current even in the absence of phase difference 
may be observed. In addition, this tangential current can flow in the opposite direction 
to the external transport current and this factor can produce the vortex-like currents. 
Finally, playing the role of magnetic Josephson phase (cj) — where q and $ are electric 
charge and magnetic flux respectively) by mis-orientation of superconducting a6-planes 
(pairing symmetry in the momentum space) is a reason for magnetic nature of pairing 
mechanism in the high T c superconductors which remains as an unknown and famous 
problem. 
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